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Abstract 

A precision measurement of the D° meson mass has been made using ~ 281 pb _1 of e + e~ 
annihilation data taken with the CLEO-c detector at the ^(3770) resonance. The exclusive decay 
D° -> K s <t> has been used to obtain M(D°) = 1864.847 ± 0.150(stat) ± 0.095(syst) MeV. This 
corresponds to M(D°D*°) = 3871.81 ± 0.36 MeV, and leads to a well-constrained determination 
of the binding energy of the proposed D°D*° molecule X(3872), as E b = 0.6 ± 0.6 MeV. 
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The D° (cu) and D ± (cd, cd) mesons form the ground states of the open charm system. 
The knowledge of their masses is important for its own sake, but a precision determination 
of the D° mass has become more important because of the recent discovery of a narrow 
state known as X(3872) [l|, @, 0, 0] • Many different theoretical models have been proposed 
(El, @, 0, 0] to explain the nature of this state, whose present average of measured masses is 
M(X) = 3871.2 ±0.5 MeV j9|. A provocative and challenging theoretical suggestion is that 
X(3872) is a loosely bound molecule of D° and D*° mesons (sj. This suggestion arises mainly 
from the closeness of M(X(3872)) to M(D°) + M(D*°) = 2M(D°) + A[M(L>*°) -M(D )] = 
2(1864.1 ± 1.0) + (142.12 ± 0.07) MeV = 3870.32 ± 2.0 MeV based on the PDG [9] average 
value of the measured D° mass, M(D°) = 1864.1 ±1.0 MeV. This gives the binding energy 
of the proposed molecule, £ b (X(3872)) = M(D°) + M(D*°) - M(X(3872)) = -0.9 ± 2.1 
MeV. Although the negative value of the binding energy would indicate that X(3872) is not 
a bound state of D° and D*°, its ±2.1 MeV error does not preclude this possibility. It is 
necessary to measure the masses of both D° and X(3872) with much improved precision 
to reach a firm conclusion. In this Letter we report on a precision measurement of the D° 
mass, and provide a more constrained value of the binding energy of X(3872) as a molecule. 

Several earlier measurements of the D° mass exist The only previous measurements 
in which sub-MeV precision was claimed are the SLAC measurements of e + e~ — *> ^(3770) — > 
D°W by the lead glass wall (LGW) [[(} and the Mark II [n| collaborations, and the CERN 
measurement by the NA32 experiment with 230 GeV 7r~ incident on a copper target (l2| . 
All three measurements determined the D° mass using D° — > K~ir + and D° — > K~7t + 7t~7t + 
(and charge conjugates) decays. In the SLAC measurements the beam constrained mass 
was determined as M 2 {D°) = E% eaia - p 2 D . The results were M(D°) = 1863.3 ± 0.9 MeV 
(LGW [13), and M(D°) = 1863.8 ± 0.5 MeV (Mark II [Tl|). The NA32 experiment reported 
M(D°) = 1864.6 ± 0.3(stat) ± l.O(syst) MeV from a simultaneous fit of the mass and lifetime 
of D° in the two decays, with the main contribution to the systematic uncertainty arising 
from magnetic field calibration. The PDG 0] lists the resulting average D° mass based on 
the measured D° masses as M(D°)avg=1864.1 ± 1.0 MeV. They also list a fitted mass, 
M(D°)fit=1864.5 ± 0.4 MeV, based on the updated results of measurements of D ± , D°, 
Df, D*^, D*°, and D*^ masses and mass differences. 

We analyze ~ 281 pb _1 of e + e~ annihilation data taken at the -?/>(3770) resonance at the 
Cornell Electron Storage Ring (CESR) with the CLEO-c detector to measure the D° mass 
using the reaction 

^(3770) -> D°W, D° -> K S (f), K s -> tt+tt", -> K + K~. (1) 

Our choice of the D° — > Ks4> decay mode is motivated by several considerations. Our 
determination of the D° mass does not depend on the precision of the determination of 
the beam energy. Since M(0) + M(Ks) = 1517 MeV is a substantial fraction of M(D°), 
the final state particles have small momenta and the uncertainty in their measurement 
makes a small contribution to the total uncertainty in M(D°). This consideration favors 
D° — > Ks4> decay over the more prolific decays D° —* Kit and D° — > Ktttctt, in which the 
decay particles have considerably larger momenta and therefore greater sensitivity to the 
measurement uncertainties. An additional advantage of the D° —* Ks<f) reaction is that 
in fitting for M(D°) the mass of Ks can be constrained to its value which is known with 
precision 0]. 

The CLEO-c detector [13] consists of a CsI(Tl) electromagnetic calorimeter, an inner 
vertex drift chamber, a central drift chamber, and a ring imaging Cherenkov (RICH) detector 
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inside a superconducting solenoid magnet providing a 1.0 Tesla magnetic field. For the 
present measurements, the important components are the drift chambers, which provide a 
coverage of 93% of 4tt for the charged particles. The final state pions and kaons from the 
decays of K$ and have momenta less than 600 MeV/c, and they are efficiently identified 
using measurements of track vertices and ionization loss (dE/dx) in the drift chambers. The 
detector response was studied using a GEANT-based Monte Carlo simulation (l4j |. 

We select D° candidates using the standard CLEO D-tagging criteria, which impose a 



very loose requirement on the beam energy constrained D° mass, as described in Ref. 15 
We select well-measured tracks by requiring that they be fully contained in the barrel region 
( | cos 1 < 0.8) of the detector, have transverse momenta > 120 MeV/c, and have specific 
ionization energy loss, dE/dx, in the drift chamber consistent with pion or kaon hypothesis 
within 3 standard deviations. For the pions from Kg decay, we make the additional require- 
ment that they originate from a common vertex displaced from the interaction point by more 
than 10 mm. We require a K$ flight distance significance of more than 3 standard devia- 
tions. We accept Kg candidates with mass in the range 497.7 ± 12.0 MeV. In addition, for 
the Ks candidates from the exclusive reaction D° — > Ks<f), we perform a mass-constrained 
(1C) kinematic fit and accept in our final sample Ks with x 2 < 20. The tt + it~ invariant 
mass distribution is shown in the upper panel of Fig. 1 with a fit to a sum of two Gaussians. 
The fit results are: M(K S ) =497. 545 ± 0.112 MeV, x 2 /d.o.f. = 0.6, and full width at half 
maximum, FWHM = 5.0 MeV. While the fit is very good, because of the limited statistics 
the resulting M(K$) does not have the precision required for testing the calibration of the 
detector. As described later, we use the large statistics data for the inclusive K s produc- 
tion, D — > Ks + X, for that purpose. The lower panel of Fig. 1 shows the K + K~ invariant 
mass distribution. The data are fitted with a Breit- Wiener of width T(0) = 4.26 MeV Q 
convoluted with the Monte Carlo determined Gaussian with FWHM = 2.8 MeV, and a 
linear background. The fit results in M(<p) = 1019.518 ± 0.243 MeV, x 2 /d.o.f. = 1.1. We 
select events containing a by requiring that M(K + K~) of the candidate kaons is within 
±15 MeV of the value M(<p) = 1019.46 MeV @. 

Figure 2 shows the invariant mass spectrum of the D° candidates constructed with Ks 
and <p as identified above. A likelihood fit of the data in the region 1840-1890 MeV was 
done with a Gaussian peak and a constant background. An excellent fit is obtained with 
the number of fitted events N(D°) = 319 ± 18, a = 2.52 ± 0.12 MeV (FWHM = 5.9 MeV), 
X 2 /d.o.f. = 0.7, and 

M(D°) = 1864.847 ± 0.150(stat) MeV. (2) 

The key to the precision measurement of the D° mass is in determining the accuracy in 
the detector calibration which can be studied by constructing M(K$) and M(<j)) from the 
measured momenta of the final state particles, and K ± . We find that M(0) is not very 
sensitive to these variations, because the K^ have very small momenta in the rest frame 
of the cf). On the other hand, M(Ks) is quite sensitive to the uncertainty in the relatively 
larger momenta of ir^ 1 in the rest frame of the Ks- The sensitivity of M(D°) is also large 
as a consequence of the sensitivity of M(K S ). We therefore conclude that M(Ks) can be 
best used to determine the accuracy of the detector calibration. As mentioned before, the 
exclusive sample of D° — > Ks4> events does not yield a statistically useful result for M(Ks)- 
It is possible to determine M(Ks) with much higher statistical precision using inclusive Ks 
production in D decays, D — > Ks + X. Inclusive Ks's were selected from each event that 
had at least one candidate D decay. The Ks mesons from the decays D° — > Ks<f) have 
momenta in the range of p{Ks) ~ 0.40 — 0.65 GeV/c. We therefore determine M(Ks) for 
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FIG. 1: Upper plot: Invariant mass of the (ir + ir~~) system for K$ decay candidates. The curve 
shows the fit with the peak shape given by the sum of two Gaussians. Lower plot: Invariant mass 
of the (K + K~) system. The curve shows the fit with a Breit-Wigner convoluted with a Gaussian 
shape and a linear background. 



this range of p(K$) in the inclusive decays. 

Figure 3 shows the M(ir + ir~) distribution for the inclusive reaction, with p(K s ) in the 
range 0.40 — 0.65 GeV/c. A fit with the peak shape given by the sum of two Gaussians and 
a linear background returns 

M{K S ) = 497.648 ± 0.007(stat) MeV. (3) 

The fit has 115,235 ± 450 events, ^jd.o.f. = 1.07, and FWHM = 4.7 MeV. 

In order to estimate the systematic error in the above determination of M(K S ), we have 
studied the variation of M(K S ) as a function of several observables associated with K s : 
p(n ± ,K s ), Pr(7r ± ), Pl^^), flight distance's), flight significance's), cos(9)(n ± , K s ), and 
7r + 7r~ opening angle. The largest variation in M(Ks) was found with respect to the variation 
in cos(#) and pt of tt + . The observed variations contribute a ±28 keV systematic uncertainty 
in our determination of M(K S ). 
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FIG. 2: Invariant mass of K$K + K system for D° — ► -ftTsc/> decay candidates. The curve shows fit 
results with a Gaussian peak shape and a constant background. 
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FIG. 3: Invariant mass of (7r + 7T~) system for K$ decay candidates from inclusive sample. The curve 
shows fit results with the peak shape given by the sum of two Gaussians, and a linear background. 

It is found that Monte Carlo events have a reconstructed output M(K S ) which differs by 
±21 keV from the input value of M(Ks). In addition, we determine systematic uncertainties 
for different peak fitting procedures: ±9 keV from variation of the peak shape, ±1 keV from 
variation of bin size from 62 keV to 250 keV, and ±8 keV from variation of fitting range 
from 15 MeV to 20 MeV. Thus, added in quadrature, the total systematic uncertainty in 
M(Ks) from the inclusive data is ±37 keV, and our final result is 

M(K S ) = 497.648 ± 0.007(stat) ± 0.037(syst) MeV. 

Since M(K s )pdg = 497.648 ± 0.022 MeV, 

M{K S ) - M(K s )pdg = 0.000 ± 0.044 MeV. 
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Systematic Error(MeV) 



Detector Calibration 


±0.090 


Monte-Carlo input / output 


±0.022 


Bin size (0.002-2 MeV) 


±0.018 


Unbinned fit 


±0.007 


Peak Shape(single/double Gaussian) 


±0.003 


Background Shape (const. /linear) 


±0.007 


Fit interval (±20 MeV) 


±0.002 


Sum in Quadrature 


±0.095 



TABLE I: Summary of systematic errors in M(D°). 



To be conservative, we consider the above maximum difference ±44 keV to be a reflection 
of the possible uncertainty in the momentum calibration of the detector, which likely arises 
from uncertainty in the magnetic field calibration and uniformity. The S-field of the CLEO- 
c detector is set by scaling a map of the -B-field such that the measured mass of J/ip — > ji + n~ 
lies at the mass of J / if) [9j . We have tried several different ways to impose ad-hoc changes in 
the measured momenta of the pions to produce a ±44 keV change in M(Kg) in the inclusive 
data. We find that when these same changes are applied to the measured momenta of all 
tt^ and K in the exclusive data, in all cases the change in M(D°) is nearly twice as large 
as the change in M(Ks). We therefore assign ±90 keV as the uncertainty in M(D°) due to 
the uncertainty in the momentum calibration of the detector. 

An independent confirmation of this conclusion is obtained by measuring the mass of 
ip(2S) via the reaction if)(2S) — > 7t + tt~J/iI), which produces tt^ with nearly the same 
momenta as 7r + 7r~ and K + K~ from the D° exclusive data. A sample of CLEO-c data 
for ip(2S) — > 7i + 7i^J/ip, J/ip — > was analyzed with the track selection and fitting 

procedure similar to those used to determine M(D°). A mass-constrained kinematic fit 
for J/tp was performed, similar to that done for the K$ in our D° decay. The fit re- 
sulted in M(i/)(2S)) = 3686.122 ± 0.021 MeV. This differs from the most precise mea- 
surement of M(ip(2S)) = 3686.111 ± 0.027 MeV by the KEDR collaboration [Tj} by 



AM(ip(2S)) = 11 ± 34 keV. Since the detector I?-field was calibrated at J/ip, this dif- 
ference can be attributed to the uncertainty in measurement of 7r + 7r~ momenta, just as in 
the case of tt + tc~ in inclusive K$- This assures us that our assignment of ±90 keV as the 
systematic uncertainty in M(D°) due to detector calibration is conservative. 

Other contributions to systematic errors in M(D°) are smaller, and are listed in Table I. 

Thus, our final result is 

M(D°) = 1864.847 ± 0.150(stat) ± 0.095(syst) MeV. (4) 

Adding the errors in quadrature, we obtain 

M(D°) = 1864.847 ± 0.178 MeV. (5) 

This is significantly more precise that the current PDG average j^]. 

Our result for M(D°) lead s to M(D°lD^) = 3871.81 ± 0.36 MeV. Thus, the binding 
energy of X(3872) as a D°D*° molecule is E b = (3871.81 ± 0.36) - (3871.2 ± 0.5) = ±0.6 ± 
0.6 MeV. This result provides a strong constraint for the theoretical predictions for the 
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decays of X(3872) if it is a D°D*° molecule [8j. The error in the binding energy is now 
dominated by the error in the X(3872) mass measurement, which will hopefully improve 
as the results from the analysis of larger luminosity data from various experiments become 
available. 
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luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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